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Abstract 

The strong coupling constants gs Q A Q7 r (Q = b and c) are studied in 
the framework of the light cone QCD sum rules using the most general 
form of the baryonic currents. The predicted coupling constants are 
used to estimate the decay widths for the Eg — ► Aqtt decays which 
are compared with the predictions of the other approaches and existing 
experimental data. 
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1 Introduction 



Recent years have witnessed advances in the heavy baryon spectroscopy, with 
the discoveries of the heavy baryons involving the b and c quarks. Since the 
spin of the baryon carries information on the spin of the heavy quark, the 
study of the heavy baryons might also lead us to study the spin effects at 
the loop level in the standard model. 

To study the meson-baryon couplings, a non-perturbative method is needed. 
Among all non-perturbative approaches, the QCD sum rules approach [T]- [3] 
has received special attention to study the properties of hadrons. In the case 
of the light baryons, this method has been successfully applied for calculation 
of the meson-baryon coupling constants. The pion-nucleon coupling constant 
has been studied in traditional three-point QCD sum rules [%]-[12j. The kaon- 
baryon coupling constants have also been calculated in the same framework 
in [13]- [IS]- The latter has also been studied in light cone QCD sum rules 
(LCQSR) in [IT]. The coupling constant for K meson-octet baryons and ir 
meson-octet baryons have also been calculated in [18] in LCQSR. 

The QCD sum rules is also applied to the study of the heavy hadron mass 
spectrum (see e.g. [19]). The masses are also studied in QCD string model 
[20] and using quark model in J2TJ[22]. In [22j, sum rules between the masses 
of the heavy baryons derived using the quark model has been analyzed and 
experimental tests of sum rules for heavy baryon masses have been discussed 
in [23]. In the present work, using the general form of the current for Sq 
and Aq baryons, we calculate the <?e q Aq7t (Q = b and c) coupling constants 
in the framework of the LCQSR approach. Having computed the coupling 
constants, we also evaluate the total decay widths for strong Eg — > Aqti 
decays and compare our results with the predictions of the relativistic three- 
quark model (RTQM) [23], light-front quark model (LFQM) [25] and existing 
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experimental data. The paper encompasses three sections: in the next sec- 
tion, we calculate the LCQSR for the coupling constant ^SqAqtt- Section III 
is devoted to the numerical analysis of the coupling constant ^EqAqtt, our 
prediction for the total decay rates and discussion. 

2 Light cone QCD sum rules for the coupling 
constant gz Q A Q n 

To calculate the coupling constant gs Q A Q 7r in LCQSR, one starts with a 
suitably chosen correlation function. In this work, the following correlation 
functions is chosen: 

n = ij d 4 xe^(n(q) | r{^ AQ (x)^ Q (0)} I 0), (1) 

where r]^ Q and ^Aq are the interpolating currents of the heavy baryons £q 
and Aq. In this correlator, the hadrons are represented by their interpolating 
quark currents. This correlation function can be calculated in two different 
ways: on the one hand, inserting complete sets of hadronic states into the 
correlation function, it can be expressed in terms of hadronic parameters such 
as the masses, residues and the coupling constants. On the other hand, it can 
be calculated in terms of quark-gluon parameters in the deep Euclidean region 
when p 2 — > — oo and (p + q) 2 — > — oo. The coupling constant is determined by 
matching these two different representations of the correlation function and 
applying double Borel transformation with respect to the momentum of both 
hadrons to suppress the contributions of the higher states and continuum. 

The derivation of the physical (or phenomenological) representation of 
the correlation function follows the same lines as in the case of light hadrons 
(see e.g. [IS])- For completeness, we repeat the derivation below. First, one 
inserts two complete sets of states between the interpolating currents in (pQ) 
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with quantum numbers of the Eg and Aq baryons. 

n (0 I Va q I A Q (p 2 )) <Sg(pi) 1 ^s Q I 0) 

n = 2 2 a q p 2 K ? £ o bi 2 — + 

Pi -m 2 A w v Pi-mi 



(2) 



where pi = p + g, P2 = P, and ... stands for the contributions of higher states 
and continuum. The vacuum to baryon matrix element of the interpolating 
currents are defined as 

(0 | T] B | B(p, s)) = X B u B (p, s), (3) 

where B = Eq or Aq, ub(p,s) is a spinor describing the baryon B and A# 
is the residue of the B baryon. The last ingredient is the matrix element 
(AQ(p2)n(q) | Eq(pi)) which can be parameterized in terms of the coupling 
constant g^ Q A Q7T as 

(A Q {p 2 )n(q) | E Q (pi)) = gx Q A Q7T u{P2)il5u{pi). 

(4) 

Using Eqs. ([2j|4]) and summing over the spin of the baryons, the following 
representation of the correlator for the phenomenological side is obtained: 

„ . ^SqAqTtAAqAsq r 

n = i— 2 r^r~2 rr i * ^ 7s ~ ms « ^ 7s 

+ {m Kq - m EQ ) ^75 + (m SQ m AQ - p 2 )-f 5 ] . (5) 

Note that, the structures ^75 and 75 have very small coefficient due to the 
fact that rn^ Q ~ m\ Q1 hence they will not yield reliable sum rules. 

To calculate the representation of the correlation function, Eq. ([1]), from 
QCD side, we need the explicit expressions of the interpolating currents for 
Eq and Aq baryons. In principal, any operator having the same quantum 
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numbers as the corresponding baryon can be used. It is well known that 
there is a continuum of choices for the heavy spin-| baryons interpolating 
currents that does not contain any derivatives. The general form of the Sq 
and Aq currents can be written as (see also 



-^=e abc {(u aT CQ b ) l5 d c + (3(u aT C l5 Q b )d c 
[(Q aT Cd b ) l5 u c + f3(Q aT C l5 d b )u c ]}, 



= -j=e abc \2[(u aT Cd b ) lb Q c + (3'{u aT C l5 d b )Q c \ + (u aT CQ b ) l5 d c 
+ (3\u aT C l5 Q b )d c + (Q aT Cd b h 5 u c + (3\Q aT C l5 d b )u c } , 



(6) 



where (3 and (3' are arbitrary parameters. For simplicity, we assume (3 = (3'. 
The (3 = —1 corresponds to the Ioffe current and C is the charge conjugation 
operator and a, b and c are color indices. 

After contracting out all quark pairs in Eq. (TjQ), the following expression 
for the correlation function in terms of the quark propagators is obtained 



n 



d*xe^(n(q) \ < y 5 S% S'™ Sf 75 



- lh S cb S' ba S a u c 75 - l/2( l5 ST S' bb S? 75 - 1^ 57 S bc 75 
+ Tr[S b Q a 'Sf} l5 Sf l5 - Tr[S b /S'f} l5 Sf^) 

+ (3 l5 S^ l5 SfS b /- l5 S§' l5 S^Sf+S^'Sf l5 S b / l5 

- S$S' ba ' l5 Sf l5 + l/2( l5 Sf l5 S'^S b / - l5 Sf l5 S' bb 'Sf 

- ST'S' bb ' l5 Sf l5 + SfS'^' l5 S b / l5 - ST'l 5 Tr[ l5 S b Q a 'S': b '} 

+ Sf l5 Tr[ l5 S b /S'f) - l5 SfTr[S b Q a ' l5 S': b '] + 75 S^' Tr[S b / lb Sf]) 



+ P 2 



^bc' 



icb' 



tlba 



cb> 



i/aa' 



ibc' 
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- Sf^S™^ - Tr[ l5 S a u b ' l5 S^']Sr + Tr[ l5 S a Q b ' 75 ST']S: 



iba' 



tab' 



0), 
(7) 



where S' = CS T C and Sq^) (q = u, d) is the full heavy (light) quark 
propagator. Note that, the Eq. (JT)) is a schematical representation for the 
full expression. To obtain the full expression from the Eq. (j7j), one should 
replace S u by u(0)u(x) to calculate the emission from the u quark, and then 
add to this the result obtained by replacing Sd by d(0)d(x). From Eq. (j7]), 
it follows that the expression of the light and heavy quark propagators are 
needed. 

The light cone expansion of the quark propagator in the external field is 
calculated in [27]. The propagator receives contributions from higher Fock 
states proportional to the condensates of the operators qGq, qGGq and qqqq, 
where G is the gluon field strength tensor. In this work, we neglect contri- 
butions with two gluons as well as four quark operators due to the fact that 
their contributions are small [28J. In this approximation, the heavy and light 
quark propagators have the following expressions: 

d 4 k 



S q( x ) = S Q ee ( x ) 



-ikx 



dv 



ft + rriQ 



[m. 



~2\2 



G» v (vx)a 



+ 



1 



S q (x) 



m. 



S{ ree { 



-i>9s 



k 2 



vx.G^ lv 



x 



x) — 
1 



du 



( gg) 

12 192 

* G, 



m (qq) 



The expression of the free light and heavy quark propagators in the x 
representation are: 



ux^G ^ v {ux)^ 1 



• (8) 



27T 2 X 4 ' 
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m 2 Q K^rngV^) _ .m% $ 
4tt 2 



4ir 2 x 2 



(9) 



where Ki are the Bessel functions. 

In order to calculate the contributions of the pion emission, the matrix 
elements (ir(q) \ qTiq | 0) are needed. Here, Tj is any member of the complete 
set of Dirac matrices {1, 75, j a , 2757a, a a p/\/2}. These matrix elements are 
determined in terms of the pion distribution amplitudes (DA's) as follows 
[291 EQ]. 
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(P)l?(a07„7sg(0)|0) = -if wPfl I due^ x [^(u) + —mlx 2 k(u 



/() 



(7r(p)|g(x)z 75 g(0)|0) 
(ir(p)\q(x)cr a ^ 5 q(0)\0} 
{7r(p)\q(x)a^ 5 g s G a/3 (vx)q(0)\0) 



2 px J 

/ due tupx (p P (u), 
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fi n (l- fil) {p a x p - p p x a ) \ due mpx (p a (u) 



PaPn 9u/3 {PvXp + Pp%i 

1 px 

1 



PaPu 9n(3 (zyC/3 + 

\ px 
\ px 



(7r{p)\q(x)^ 5 g 8 G a0 (vx)q(O)\O) 



= P»{p a x p -ppx a )—Uml [ Vae^ +va ^ x A\\ 
px J 



+ 



Pf3 ( 9^a (P^X a +p a X^ 
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- Pa \9nP ~ +Pf3 X V^J f* m l 

(7i(p)\q(x) lfl tg s G a p(vx)q(0)\0) = Plt (p a x p - p P x a )— f v m* ( Vae l ^ +va ^ x V\\ 

px J 

+ P/3 \ 9na - ~(P» X <* + Pa^Sj 

x / Vae i ( a * +va « )px V±(a i ), 



on 



- Pa [ 9^ - —(P^ X P + PP X n) 



where fx n = f n 



m u +m d 

T7W 



and the functions <p n (u) 



" r m u +m ( j' 

<Pp(u), ipo-(u), T(cKi), Al(atj), Vj_(a;i) and V||(oti) are functions of 

definite twist and their expressions will be given in the numerical analysis 
section. The measure Va is defined as 



T>a = / docq I da q I da g 5(l — a q — a q — a g ). (11) 
Jo Jo Jo 

Note that, in the approximation of this work where we neglect the light quark 
masses, m 2 = 0, fi n = 0, = —(uu)/f w = —(dd)/f n . 

Using the expressions of the light and heavy full propagators and the pion 
DA's, the correlation function Eq. ([1]) can be calculated in terms of QCD 
parameters. Separating the coefficient of the structure p 1 ^75 in both repre- 



sentations, and equating them, sum rules for the coupling constant 
is obtained. The contribution of the higher states is subtracted using quark 
hadron duality, and in order to further suppress their contribution, Borel 
transformation with respect to the variables p\ = p 2 and p\ = (p + q) 2 is 
applied. Here, we should mention that we have also studied the other struc- 
ture in Eq. (j3J), i.e., ^75 but its result for coupling constant is not stable and 
only the ^75 structure leads to reliable prediction on the coupling constant 



(10) 
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S'SqAqtt- 

The sum rules for the coupling constant is obtained as 



where the function II is 



IT, 



n 



so 



— s ^ 

e* 5 p(s)ds + e~S 2 ~r, 



(12) 



(13) 



with 

P(s) 



< dd > + < uu >) 



12^6 



m Q (/3-l) -6 



2(^20 - MmQ^i-Ul + 2/3) + Ce(l + P)] 



^io"Iq^[3C B (1 + /?) - 4Ce] - m Q ^(3C 5 (l + /3) - 4C 6 )M— ) 



+ 6/X(l + /3)[2^io - ^20 + ^3i + 2/n(-!i)^(Mo) 



+ 2(^20 - i>3l)mQlJLn(l + 2/3)<f a (u Q ) 



(14) 



and 

r = 



192^ 



'< dd > + < uu >) 



^[-11 - 17/3 + (7 + /?)](/? - l)w(ito) 



rriQii^p 2 + p + l)(p a (u ) 



3M 2 

9Um Q (P ~ - ^(3/3 2 + 2/3 + 3)^ CT (n ) 

— -=(< dd> + <uu >)mQ/i 7r (/3 2 + (3 + l)^ CT (tt ). 
6v6 



+ 



The other functions entering Eqs. f 1 1 4ti 1 5 [) are given as 



Von / dvfj(ai)5(a q + va g - u ), 
o 



(15) 
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■> = I va -I! 



dvgj(ai)S'(a q + va g — ^o), 

= (g — WlQ ) 

(16) 

and = V\\(ai), / 2 («i) = vV||(o:i), h{ai) = V±(«i), / 4 («i) = vV±(ai), 

9i{oti) = T(oii) and (?2(«;) = vT(ai) are the pion distribution amplitudes. 
Note that, in the above equations, the Borel parameter M 2 is defined as 
M 2 = m^+m 11 anc ^ u ° = m^+m^ • Since the mass of the initial and final 
baryons are close to each other, we can set M\ = M| = 2M 2 and u — |. The 
contributions of the terms ~< G 2 > are also calculated, but their numerical 
values are very small and therefore for customary in the expressions these 
terms are omitted. 

For calculation of the coupling constants of the considered baryons, their 
residues, \t, q (a q ) are needed. Their expressions are obtained as: 

~ X k^) e ~ mlQ{AQ)/M2 = r ^Pi( 2 )(^ + e^r 1(2) , (17) 

Jm Q 

with 

Pl (s) = (<dd> + <uu >) ^ 64?r2 1 4^(^00 ~ 13^02 - fyu) 
+ 3m Q (2^i - - ^12 + 2^>2i) 

4 

+ ^h^ 5 + ^ 2 + M 12 ^ - 6 ^o + 2^0 - 4^41 + ^42 - 12/n(^)], 
2U487T 4 vtlq 

(18) 

p 2 (s) = (<dd> + <uu >)^— f|- 



rji 



- 6(1 + 0)ih,] + (1 + 5/?)m Q (2<h„ - Vn - V12 + 2V'2i) 



fnh s 
+ ^t#t[5 + m + 5/3) 12^o - 6^20 + 2^30 - #41 + ~ 12Zn — L 



2048vH 



(19) 



r 2 



(/?-l) : 
24 

(/?-i) 
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■ < dd >< uu > 
< dd >< uu > 



2M 4 4M 2 



m 2 Q ml 



2M 4 



(13 + 11/?) 



+ £^(25 + 23/3) -(13 + 11/3) I. (20) 

3 Numerical analysis 

This section is devoted to the numerical analysis for the coupling constant 
^SqAqtt and calculation of the total decay width for Sq — > Aq7t. The 
input parameters used in the analysis of the sum rules are (uu)(l GeV) = 
(dd)(lGeV) = -(0.243) 3 Gel/ 3 , 8(1 GeV) = 0.8(mx)(l GeV), m b = 4.7 GeV, 
m c = 1.23 GeV, m Et = 5.805 GeV, m Sc = 2.439 GeV, m Ab = 5.622 GeV, 
m Ac = 2.297 GeV, and mg(l GeV) = (0.8 ± 0.2) GeV 2 [3XJ. From the sum 
rules for coupling constant, it is clear that the 7r-meson wave functions are 
needed. These wave functions are given as 



u 



T{0Li 



) P {U) 



+ 



4>a(u) 

Vntau) 



6uu (l + alC x (2u - 1) + a£C 2 2 (2u - : 
360r] 3 a q a q a 2 g ^1 + w 3 ^(7a g - 3) J , 

1+ f307 ?3 -~^C|(2u-l) 

-3ife«* " ^ " TsA-^) C 4 (2u - 1), 

20 /4 io //J y 

(1 7 3 \ 3 

5^3 - 2^3 - — nl - -nla% j Cj(2u - V 

I20a g aqa g (v 00 + v w (3a g - 1)) , 
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A\\(<Xi) 
V ± (a z ) 

A±(ai) 
B{u) 

gM 

A(u) 



= 120a q a q a g (0 + a w (a q — a q )) , 



30a g (a q — a q ) 
g w (u) - ff (u), 



^oo(l - oig) + h 01 (a g (l - a g ) - 6a q a q ) + h w (a g (l - a g ) - -(a? + a 2 )) 

1 



ft-oo + h ia g + -h 10 (5a g - 3) 



= g C~ 2 (2u - 1) + g 2 Ci(2u - 1) + <? 4 C| (2m - 1), 



6mm 



16 24 „ on 20 
- + -4 + 20,3 + Y . 4+ - 



' - ^3^3 - ^4 ) C|(2m - 1) 



15 16 27 



27 



+ 



-^;-T^ 3 )cj(2„-i) 



+ (~Y a 2 + 21774^4^ [2m 3 (10 - 15m + 6m 2 ) Iiim 
+ 2m 3 (10 - 15m + 6m 2 ) Iiim + mm(2 + 13mm)] , 

where C*(x) are the Gegenbauer polynomials, 

1 



^00 

VlO 

hw 

go 
g2 

g± 



21 



j?74, 

9 



21 



4 7?4W4 ~ 20 a2 ' 



A VaWa + 20° 2 ' 



7 

4 
7 

4 

1, 

18 „ 20 
1 + + 60 ^3 + y774, 

9 

~28 a2 _ 6r?3W3 ' 



(22) 



The constants appearing in the wave functions are calculated at the renor- 
malization scale /j, — 1 GeV 2 and they are given as a\ = 0, a\ = 0.44, 
i] 3 = 0.015, 774 = 10, m> 3 = —3 and m> 4 = 0.2. 



(21) 
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The sum rules for the coupling constant also contains three auxiliary 
parameters: Borel mass parameter M 2 , continuum threshold s and general 
parameter (3 enters the expressions of the interpolating currents. In principal, 
M 2 and (5 are completely arbitrary and hence the coupling constant, which 
is a physical observable, should be independent of their exact values. In 
practice, though, due to the approximations made in the calculations, there 
is a residual dependence of the predictions on these unphysical parameters. 
Hence, a range for these parameter should be found where the predictions 
are practically insensitive to variations of these parameters. To find the 
working region for M 2 , we proceed as follows. The upper bound is obtained 
requiring that the contribution of the higher states and continuum should be 
less than that of the ground state. The lower bound of M 2 is determined 
from condition that the highest power of 1/M 2 be less than say 30°/ of 
the highest power of M 2 . These two conditions are both satisfied in the 
region 15 GeV 2 < M 2 < 30 GeV 2 and 4 GeV 2 < M 2 < 10 GeV 2 for baryons 
containing b and c-quark, respectively. The third parameter, s has a physical 
meaning, and it should have a value near the first excited state. The value 
of the continuum threshold is calculated from the two-point sum rules. We 
choose the interval s = (6.0 2 - 6.2 2 ) GeV 2 and s = (2.5 2 - 2.7 2 ) GeV 2 for 
baryons containing the b and c quark, respectively. 

In Figs. 1 and 2, we present the dependence of the coupling constants 
9s b A b w an d (fe c A c7r , at fixed values of the continuum threshold s and the 
general parameter (5. From these figures, we see a good stability for coupling 
constants <?e 6 a 6 7t and gT, c A c n with respect to the Borel mass square M 2 in the 
working region. The next step is to determine the working region for auxiliary 
parameter (3. For this aim, in Figs. 3 and 4, we depict the dependence of 
the coupling constants gs b A b n and gs c A c n on cosO where tan6 = j3, at two 
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fixed values of M 2 . From these Figures, we see that the best stability for the 
coupling constants gs b A b -n and (7e c a c tt is in the region —0.5 < cos9 < 0.2. 
Our final results on coupling constants g-£ b A b7T and gs c A C 7r are: 

9n b A b n = 23.5 ±4.9, 
gx c A C Tr = 10.8 ±2.2. 

(23) 

The quoted errors are due to the uncertainties in the input parameters as 
well as variation of the Borel parameter M 2 , continuum threshold So and 
general parameter (3. 

Having computed the coupling constant g-£ Q A Q n, the next step is to cal- 
culate the total decay width for — > A^n and S c — > A c tt decays. From 
Eq. (jl]) the transition amplitude is M = ^SqAqtt^TsM and the differential 
decay width is found in terms of the coupling constant as: 

12 



\9xqA 



(24) 



where \lf\ = (ra| Q — m A Q )/^ m ^Q- The numerical values of the decay rates 
are given in Table 1. In order to compare with the predictions of other 
methods, in the same table, we present the predictions of the relativis- 
ts three-quark model (RTQM) [21], light-front quark model (LFQM) [25] 
and existing experimental data [32]. This table depicts a good consistency 
among the methods and the experimental data in order of magnitudes for 
charm case. Note that, due to the isospin symmetry the decays of differ- 
ent charge — > A+7r + '°'~(£jj~' '~ — > Af } 7i +, °'^) have the same decay 
widths. Experimentally, only the widths for — > A+7r +, ~ are measured 
and the value in the table is their average. Only the upper bound for the 
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— > A+7T is known and it is consistent with other decay modes. Our 
prediction for the decay rate of the bottom case can be tested in the future 
experiments. 





T(S C — > A c tt) 


— > A b n) 


Present work 


2.16 ±0.85 


3.93 ± 1.5 


RTQM [24J 


3.63 ±0.27 




LFQM [25J 


1.555 ±0.165 




Exp. [32J 


2.21 ±0.40 





Table 1: Results for the decay rates of Eg — ► Aqir in different approaches 
in MeV. 

In summary, we calculated the gy, b \ b -K and g-£ c A c7T coupling constants in 
the light cone QCD sum rules approach. Using these coupling constants, we 
also evaluated the total decay width for the strong £& — > A&7T and S c — > 
A c 7r decays and compared with the predictions of the other approaches and 
existing experimental data. 
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Figure 1: The dependence of the gv. b K b -K on the Borel parameter M 2 at fixed 
value of the continuum threshold s n = 6.0 2 . 
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Figure 2: The same as Fig. 1, but for gs c A C 7r and fixed value of the continuum 
threshold s = 2.5 2 . 
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Figure 3: The dependence of |(fe b A b7I -| on cos9 at fixed value of the continuum 
threshold s = 6.0 2 . 




Figure 4: The same as Fig. 3, but for |<?e c a c tt| and fixed value of the contin- 
uum threshold s = 2.5 2 . 
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